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Abstract 



We propose a technique for the design and analysis of decentrahzed adaptation algo- 
rithms in interconnected dynamical systems. Our technique does not require Lyapunov 
stability of the target dynamics and allows nonlinearly parameterized uncertainties. We 
show that for the considered class of systems, conditions for reaching the control goals 
I can be formulated in terms of the nonlinear L2-gains of target dynamics of each intercon- 

' nected subsystem. Equations for decentralized controllers and corresponding adaptation 

algorithms are also explicitly provided. 

Keywords: nonlinear parametrization; unstable, non-equilibrium dynamics; decentral- 
ized adaptive control; monotone functions 

> ■ 

(N Notation 
o . 

\0 . According to the standard convention, M defines the field of real numbers and ]R>c = {x G 
M.\x > c}, ]R+ = ]R>o; symbol R"' stands for a linear space £(R) over the field of reals with 
^ dim{£(R)} = n; ||x|| denotes the Euclidian norm of x G R"; denotes the space of functions 
^ ■ that are at least k times differentiable; /C denotes the class of all strictly increasing functions 
K : R_,_ such that k.{0) = 0. By Lp[tQ,T], where T > 0, p > 1 we denote the space of 

^ ■ all functions f : R_|_ R" such that ||f ||p,[f,),T] = ( /g l|f < oo; 11^ llp,[to,T] denotes the 

' Lp[to, T]-norm of f{t). By L^l^lto, T] we denote the space of all functions f : R+ R" such that 
||f||oo,[io,T] = esssup{||f(t)||,t G [to,T]} < oo, and ||f ||oo,[to,T] stands for the L'^[tQ,T] norm of 

A function f(x) : R'^ R"^ is said to be locally bounded if for any ||x|| < 5 there exists 
a constant D{6) > such that the following inequality holds: ||f(x)|| < D{6). Let F be an 
n X n square matrix, then F > denotes a positive definite (symmetric) matrix, and F~^ is the 
inverse of F. By F > we denote a positive semi-definite matrix, ||x||p to denotes the quadratic 
form: X"^Fx, x G R". The notation | ■ | stands for the modulus of a scalar. The solution of a 
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system of differential equations x = f(x, t, 0,u), x(to) = xq, u : M+ —>■ R"", e eW^ foTt> to 
will be denoted as x(t, xq, to, ^, u), or simply as x(t) if it is clear from the context what the 
values of xq, are and how the function u{t) is defined. 

Let u : M" X R'^ X R+ ^ M"^ be a function of state x, parameters 6, and time t. Let in 
addition both x and be functions of t. Then in case the arguments of u are clearly defined 
by the context, we will simply write u(t) instead of u(x(t), ^(t), t). 

The (forward complete) system x = f (x, t, 6, u(t)), is said to have an L^[to, T] i— > L^[to, T], 
gain (T > t^, p,q E ]R>i U oo) with respect to its input u{t) if and only if x(t, Xq, to, 0, u(t)) G 
L^[to,T] for any u(t) G i^^Tf^o,^] and there exists a function 7,,^ : x x ^ R+ 
such that the following inequality holds: ||x(t)||q^[tQ^r] < 7q^p(xo, 0, ||u(t)||pjtQ^T])- The function 
7(},p(xo, 0, \\'^(t)\\p^[to,T]) is assumed to be non- decreasing in ||u(t)||pjt(, ^j, and locally bounded in 
its arguments. 

For notational convenience when dealing with vector fields and partial derivatives we will 
use the following extended notion of the Lie derivative of a function. Let x G R" and assume x 
can be partitioned as follows x = xi © X2, where xi G R^, xi = (a^n, . . . , Xig)'^, X2 G R^, X2 = 
{x2i, . . . , X2pY' 1 q + p = n, and © denotes the concatenation of two vectors. Define f : R" — R" 
such that f (x) = fi(x) © f2(x), where fi : R" ^ R^ fi(-) = (/n(-), • • • , fiq{-)f, h ■ ^ R^ 
f2(-) = (/2i(-)5 • • • ) /2p('))"^- Then Lf-(x)'?/'(x, t), i G {1,2} denotes the Lie derivative of the 
function ^/j{'K,t) with respect to the vector field fj(x, 0): Lf-(x)'?/'(x, t) = X]^'™^' 

1 Introduction 

We consider the problem how to control the behavior of complex dynamical systems com- 
posed of interconnected lower- dimensional subsystems. Centralized control of these systems 
is practically inefficient because of high demands for computational power, measurements and 
prohibitive communication cost. On the other hand, standard decentralized solutions often face 
severe limitations due to the deficiency of information about the interconnected subsystems. 
In addition, the nature of their their interconnections may vary depending on conditions in 
the environment. In order to address these problems in their most general setup, decentralized 
adaptive control is needed. 

Currently there is a large literature on decentralized adaptive control which contains suc- 
cessful solutions to problems of adaptive stabilization [HI El, tracking [71 |HJ [T71 HH], and out- 
put regulation [HI EBj of linear and nonlinear systems. In most of these cases the problem 
of decentralized control is solved within the conventional framework of adaptive stabiliza- 
tion/tracking/regulation by a family of linearly parameterized controllers. While these results 
may be successfully implemented in a large variety of technical and artificial systems, there 
is room for further improvements. In particular, when the target dynamics of the systems is 
not stable in the Lyapunov sense but intermittent, meta-stable, or multi-stable [H UHl UHl or 
when the uncertainties are nonlinearly parameterized (3 El IH and no domination of the 
uncertainties by feedback is allowed. 

In the present article we address these issues at once for a class of nonlinear dynamical sys- 
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terns. Our contribution is that we provide conditions ensuring forward-completeness, bound- 
edness and asymptotic reaching of the goal for a pair of interconnected systems with uncertain 
coupling and parameters. Our method does not require availability of a Lyapunov function for 
the desired motions in each subsystem, nor linear parametrization of the controllers. Our re- 
sults can straightforwardly be extended to interconnection of arbitrary many (but still, a finite 
number of) subsystems. Explicit equations for corresponding decentralized adaptive controllers 
are also provided. 

The paper is organized as follows. In Section 2 we provide a formal statement of the problem. 
Section 3 contains necessary preliminaries and auxiliary results. In Section 4 we present the 
main results of our current contribution, and in Section 5 we provide concluding remarks to 
our approach. 

2 Problem Formulation 

Let us consider two interconnected systems and Syi 

S^: ± = i{x,0^) +'jy{y,t) + g{x)u^ (1) 
'Sy-- y = q(y,0y) +7:r(x,)f:) + z(y)Mj^ (2) 

where x G M""", y G M"^ are the state vectors of systems and Sy, vectors Or^ G M""^ , Oy G M^^f 
are unknown parameters, functions f : M"^ x R"»^ — » M"^, q : W'y x W^^v R"-y , g : M"^ 
M"^, z : W"-y are continuous and locally bounded. Functions 7y : W^y x IR+ ^ ]R„, 

'■ IR""" X MJ^y , stand for nonlinear, non-stationary and, in general, unknown couplings 

between systems Sx and Sy, and G M, Uy G M are the control inputs. 

In the present paper we are interested in the following problem 

Problem 1 Let ipx : M"^ x IR+ M, -(/^j^ : MJ^y x IR+ — > R be the goal functions for systems 
Sx, Sy respectively. In the other words, for some values Ex G IR+, Sy G M+ and time instant 
t* G M+, inequalities 

< ^x, \\4jy{y{t),t)\\ 

< (3) 

specify the desired state of interconnection (H}, Derive functions Ux(x, t), Uy{y, t) such that 
for all Ox G M"»- , Oy G R'^'y 

1) interconnection (0), ((21) is forward-complete; 

2) the trajectories x(t), y(t) are bounded; 

3) for given values of Ex, Sy, some t* G exists such that inequalities © are satisfied or, 
possibly, both functions iljx{'^{t),t), ipy{y{t),t) converge to zero as t ^ oo. 

Function Ux{-) should not depend explicitly on y and, symmetrically, function Uy{-) should 
not depend explicitly on x. The general structure of the desired configuration of the control 
scheme is provided in Figure 1. 
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Figure 1: General structure of interconnection 



In the next sections we provide sufficient conditions, ensuring solvability of Problem ^ and 
we also explicitly derive functions ^^^(x, t) and Uy{y,t) which satisfy requirements 1) - 3) of 
Problem ^ We start with the introduction of a new class of adaptive control schemes and 
continue by providing the input-output characterizations of the controlled systems. These 
results are given in Section 01 Then, using these characterizations, in Section 0] we provide the 
main results of our study. 

3 Assumptions and properties of the decoupled systems 

Let the following system be given: 



(4) 



Xi =fl(x) + gl(x)M, 
±2 =f2(x, 0) +g2(x)M, 

where 

xi = (xii, . . . , xiqf e R"; X2 = (X2i, . . . , X2pf G MT; 

X = {Xu, Xiq, X21, . . . , X2p)^ G 

6 ^ Qg E M.'^ is a vector of unknown parameters, and Qg is a closed bounded subset of M.'^; 
M G M is the control input, and functions fi : M" ^ fa : M" x ^ W, gi : R" ^ R^, 
g2 : R" MP are continuous and locally bounded. The vector x G R" is the state vector, and 
vectors Xi, X2 are referred to as uncertainty-independent and uncertainty- dependent partition 
of X, respectively. For the sake of compactness we will also use the following description of Q: 

i = f(x,0) + g(xK (5) 

where 

g(x) = (^ll(x), . . . , ^ig(x), ^2l(x), . . . , ^2p(x))'^, 

f (x) = (/n(x), . . . , /i,(x), /2i(x, 0), . . . , /2,(x, 0)f. 
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As a measure of closeness of trajectories x(t) to the desired state we introduce the error or 
goal function ip : MJ^ x M+ — R, G C^. We suppose also that for the chosen function ip{^,t) 
satisfies the following: 

Assumption 1 (Target operator) For the given function il){^,t) G the following property 
holds: 

||x(t)||oo,[to,T] < 7 (XO,^, ||^(x(t),t)||oo,[to,T]) (6) 

where 7 (xo,0, t)||oo,[to,T]) is a locally hounded and non-negative function of its argu- 

ments. 

Assumption^ can be interpreted as a sort of unboundedness observability propeTtj |10j of system 
(jH) with respect to the "output" function ?/'(x, t). It can also be viewed as a bounded input 
- bounded state assumption for system (jH) along the constraint '?/'(x(t, xq, to, ^, w(x()f:), t)), t) = 
v{t), where the signal v{t) serves as a new input. If, however, boundedness of the state is 
not explicitly required (i.e. it is guaranteed by additional control or follows from the physical 
properties of the system itself). Assumption [T] can be removed from the statements of our 
results. 

Let us specify a class of control inputs u which can ensure boundedness of x(t, xq, to, "w) 
for every 6 E Qg and xq G M". According to (jB)), boundedness of x(t, xq, to, ^5 w) is ensured 
if we find a control input u such that '?/^(x(t), t) G Ll^[to,oc]. For this objective consider the 
djTiamics of system ^ with respect to ^/'(x, t): 

dip(-x,t) 

i' = I'f(x,0)^/'(X, t) + Lg(x)Z/'(x, t)u H C^) 

Assuming that the inverse (Lg(x)'0(x, t)) ^ exists everywhere, we may choose the control input 
u in the following class of functions: 

1 / dipha t)\ 

^' = (,-^^(-^)^("' - ^' - ) (8) 

where a; G fi^ C is a vector of known parameters of the function ip{ip,uj,t). Denoting 

Lf(x,6»)^(x, t) = f{'x,6,t) and taking into account (jH} we may rewrite equation in the 
following manner: 

^ = /(x, e, t) - /(x, e, t) - <^(^, a;, t) (9) 

For the purpose of the present article, instead of Q it is worthwhile to consider the extended 
equation: 

^ = /(x, 0, t) - /(x, e, t) - ip{,p, cv, t) + e{t), (10) 

where, if not stated overwise, the function e : R_|_ ^ M, e G Ll[tQ, oo]r\C^. One of the immediate 
advantages of equation (fTTH) in comparison with is that it allows us to take the presence of 
coupling between interconnected systems into consideration. 
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Let us now specify the desired properties of the function '{){;ip, w, t) in (jHJ, ((TUJ. The majority 
of known algorithms for parameter estimation and adaptive control ^21 ^1 ^1 assume 
global (Lyapunov) stability of system for 6 = 6. In our study, however, we refrain from 
this standard, restrictive requirement. Instead we propose that finite energy of the signal 
f{'x{t),6,t) — f{-x(t),6{t),t), defined for example by its Ll[to,oo] norm with respect to the 
variable t, results in finite deviation from the target set given by the equality ip{yi,t) = 0. 
Formally this requirement is introduced in Assumption |21 

Assumption 2 (Target dynamics operator) Consider the following system: 

^ = -^{^,u;,t) + at), (11) 

where ( : R_|_ R and (f{ip,LO,t) is defined in mU\} . Then for every uj & system ill]) has 
Ll[tQ,oo] \—>- L^[to,oo] gain with respect to input C{t)- Iit- other words, there exists a function 
7oo,2 such that 

\\m\\oo,[to,T] < 7oo,2(^o,^, ||CWi|2,[to,T]), V at) G Ll[to,T] (12) 

In contrast to conventional approaches. Assumption |21 does not require global asymptotic sta- 
bility of the origin of the unperturbed (i.e for ((t) = 0) system (fTTj) . When the stability of the 
target dynamics = —ip{ip, uj, t) is known a-priori, one of the benefits of Assumption El is that 
there is no need to know a particular Lyapunov function of the unperturbed system. 

So far we have introduced basic assumptions on system (0} and the class of feedback con- 
sidered in this article. Let us now specify the class of functions /(x, 6,t) in |TTH) . Since general 
parametrization of function f{-x,6,t) is methodologically difficult to deal with, but solutions 
provided for nonlinearities with convenient linear re-parametrization often yield physically im- 
plausible models and large number of unknown parameters, we have opted for a new class of 
parameterizations. As a candidate for such a parametrization we suggest nonlinear functions 
that satisfy the following assumption: 

Assumption 3 (Monotonicity and Growth Rate in Parameters) For the given function 
f{-x,6,t) in nj^) there exists function q;(x, t) : x M_|_ —>■ M'', a(x, t) G and positive con- 
stant D > such that 

(/(x, 6, t) - /(x, 6, t))(a(x, tfi6 - 6)) > (13) 

|/(x, 6, t) - /(x, 6, t) I < Z}|a(x, tf{6 - 6) \ (14) 

This set of conditions naturally extends from systems that are linear in parameters to those 
with nonlinear parametrization. Examples and models of physical and artificial systems which 
satisfy Assumption El (at least for bounded 6,6 & Qq) can be found in the following references 
[2113111111 mi- Assumption El bounds the growth rate of the difference |/(x, 6, t) — /(x, 6, t) \ 
by the functional D\cx{'x,t)'^ {6 — 6)\. In addition, it might also be useful to have an estimate 
of |/(x, 6, t) — /(x, 6, t)\ from below, as specified in Assumption El 
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Assumption 4 For the given function f{'x,6,t) in ilU^) and function ct{'x,t), satisfying As- 
sumption [3 there exists a positive constant Di > such that 

|/(x, e, t) - /(x, e,t)\> A |a(x, tf{e - e) \ (is) 

In problems of adaptation, parameter and optimization estimation, effectiveness of the algo- 
rithms often depends on how "good" the nonlinearity /(x, is, and how predictable is the 
system's behavior. As a measure of goodness and predictability usually the substitutes as 
smoothness and boundedness are considered. In our study, we distinguish several of such spe- 
cific properties of the functions f{'x,0,t) and {p{ilj,u),t). These properties are provided below. 

H 1 The function /(x, 6,t) is locally bounded with respect to x, 6 uniformly in t. 

H 2 The function /(x, G C^, and df{'x,6,t)/dt is locally bounded with respect to x, 6 
uniformly in t. 

H 3 The function (f{ilj,u:,t) is locally bounded in ip, uj uniformly in t. 

Let us show that under an additional structural requirement, which relates properties of the 
function q;(x, t) and vector-field f(x, 0) = fi(x, 0) © f2(x, 0) in (jlj), there exist adaptive 
algorithms ensuring that the following desired property holds: 

x(t) gL^ [to, oo]; /(x(t),0,t)-/(x,^(t),t) GLi[to,oo] (16) 

Consider the following adaptation algorithms: 

^(x, t) = r(^p(x, t) + Oiit)); r G M'^^'^, r > 

6»p(x, t) = ^(x, t)Q:(x, t) - ^(x, t) (17) 
Oi = <^(^(x, t), u), t)Q:(x, t) + 7^(x, e, u(x, 0, t),t), 

where the function 7^(x, 6, n(x, 0, t), t) : M" x M'^ x M x M+ ^ M"^ in (dH) is given as follows: 
n{x,u{yi,e,t),t) = 9^(x,t)/9t - ^/'(x,t)(9Q:(x,t)/(9t + Lf,CK(x,t)) 

(18) 

+ Lf,^(x,t) - (V^(x,t)Lg,a(x,t) - Lg,^(x,t))M(x,6»,t) 
and function ^(x,t) : M" x M+ ^ M^, ^(x,t) G satisfies Assumption El 
Assumption 5 There exists a function \l/(x, t) such that 

Additional restrictions imposed by this assumption will be discussed in some details after we 
summarize the properties of system (jH), (jHI), (fTTj) . (fTH|) in the following theorem. 
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Theorem 1 (Properties of the decoupled systems) Let system / flD)) . (fTTp , / f7^) &e 

g'ifen and Assumptions\^ &e satisfied. Then the following properties hold 

PI ) Let for the given initial conditions x(to), Oiito) and parameters vector 6, interval [to, T*] 
be the (maximal) time-interval of existence of solutions of the closed loop system / filTj) . p7[ ), 
(C^j. Then 

||/(x(t),0,t)-/(x(t),^(t),t))||2,[i„T*] <I^/(^,to,r,||e(t)||2,[to,Ti); (20) 

D^(0,to,r,ii5(t)ii2,[to,T.]) = l^j\\0 - mnu^ + §-Mt)h,[to,T'] 

\\e - 0{t)\\l-, < \\e{to) - e\\l-, + ^ "-^-^"^ 



2^lFl''^ll2,[io,T*] 

In addition, if AssumptionsU\ and{^ are satisfied then 
P2) ^{^{t),t) G Ll[to,oo], x(t) G L^[to,oo] and 

llV^(x(t),t)|| oo,[to,oo] ^ Too, 2 

(V^(xo,to),a;,P) (21) 

V = Dj{e,to,T, ||£(t)||2,[to,oo]) + lk(t)||2,[te,oo] 

PS'y) if properties iOl /ioW, anc? system ill]) has Ll[tQ, oo] i— > -^Vp[to, C)o], p > 1 gain with 
respect to input ({t) and output ip then 

e{t) G L^[to,oo] n L^[to,oo] ^ lim ^(x(t),t) = (22) 

//, m addition, property I^holds, and the functions a(x, t), dip{^,t) / dt are locally bounded 
with respect to x uniformly in t, then 
P4 ) the following holds 

lim /(x(t), e, t) - /(x(t), 0{t),t) = (23) 

t— +00 

The proof of Theorem ^ and subsequent results are given in Section 6. 

Let us briefly comment on AssumptionEl Let Q:(x, t) G C^, Q:(x, t) = col(ai(x, t), . . . , arf(x, t)), 
then necessary and suflicient conditions for existence of the function \l/(x, t) follow from the 
Poincare lemma: 

« ^,i.,tf-^]^(^(,.(.,tf-^)Y (24) 



9X2 \ 9X2 / V9X2 \ 9X2 

This relation, in the form of conditions of existence of the solutions for function \E'(x, t) in 
takes into account structural properties of system (jlj), (fTn|) . Indeed, consider partial derivatives 
9aj(x, t)/9x2, dip(x,t)/d:K2 with respect to the vector X2 = (x2i, . . . ,a;2p)"^. Let 

^^=(0 ... * ... 0),^^=( ... * ... ) (25) 

9X2 9X2 

where the symbol * denotes a function of x and t. Then condition ()25|) guarantees that equality 
(and, subsequently. Assumption holds. In case 9a(xi ©X2,t)/9x2 = 0, Assumption 
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holds for arbitrary ilj{'x,t) G C^. If ip{'x,t), q;(x, t) depend on a single component of X2, for 
instance X2k, k E {0, . . . ,p}, then conditions hold and the function \I'(x, t) can be derived 
explicitly by integration 

^(x,t) = /"^(x,t)^^rfx2fc (26) 

J OX2k 

In all other cases, existence of the required function \l/(x, t) follows from ()24|) . 

In the general case, when dim{x2} > 1, the problems of finding a function \l/(x, t) satisfying 
condition (fTlH) can be avoided (or converted into one with an already known solutions such as 
()24|1 . (j2ni)) by the embedding technique proposed in [21]] • The main idea of the method is to 
introduce an auxiliary system that is forward-complete with respect to input x(t) 



(27) 



^ = f^(x,tt), iew 

such that 

||/(x(t),0,t) - /(xi(t) ©h5(t) ©x'2(t),0,t)||2,[io,T] < Q e R+ (28) 
for all T > to, and dim{h^} + dim {X2} = p. Then can be rewritten as follows: 

i) = /(xi © © x^, 6>, t) - /(xi © © x^, 0, t) - (p{ip, LJ, t) + £^{t), (29) 

where e^{t) G Ll[tQ, 00], and dim{x2} = p — h < p. In principle, the dimension of Xg could be 
reduced to 1 or 0. As soon as this is ensured. Assumption El will be satisfied and the results 
of Theorem ^ follow. Sufficient conditions ensuring the existence of such an embedding in the 
general case are provided in [2^. For systems in which the parametric uncertainty can be 
reduced to vector fields with low-triangular structure the embedding is given in [21]. 

4 Main Results 

Without loss of generality let us rewrite interconnection (QJ), Q as follows : 

Xi = fi(x) + gi(x)M^ 

(30) 

X2 = f2(x, 0^) + 7y(y, t) + g2(x)M3, 



yi = qi(y) + zi(y)% ^^^^ 
y2 = q2(y, Oy) + 7^(x, t) + Z2iy)uy 

Let us now consider the following control functions 



u 



^(X, 0^, U}^, t) = (Lg(x)V^x(x, t)) ^ (^-Lf(^_^^)?/'^(X, t) - (fi^ilpcc, i^x, t) 

R.rxE,^R (32) 
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(34) 



Uy{y, 0y, UJy, t) = (L^(y)^y(y, t)) ^ (^-i^q(y,0^)^y(y, t) - (fyilpy, UJy, t) 

dipJy,t)\ (33) 

— dt ) ' ■ ^ ^ ^+ ^ ^ 

These functions transform the original equations (jHUj) . (jHT|) into the following form 

4'x = -^xitpx, ^x, t) + /x(x, 9^, t) - /^.(x, e^, t) + /iy(x, y, t) 
i^y = -^yi-ipx, ^y, t) + fy{y, Oy, t) - fy{y, Oy, t) + K{x, y, t), 

where 

^x(x,y,t) = L^^(y,t)?/^^(x,t), /iy(x,y,t) = L^4x,t)V'?;(y, ^) 
/^(x,0^,t) = Lf(x,0,)^/':,(x,t), fy{yi,Oy,t) = Lq(y,0^)^/;y(y,t) 
Consider the following adaptation algorithms 

^,(x,t) = r,(6/p,,(x,t) + ^/,..(t)); G M'^^^ r, > o 

^P,^.(x, t) = Vx(x, t)a^.(x, t) - ^^.(x, t) (35) 
= V2x(V'x(x, t),u;.^, t)cx.^{x, t) + 7^^.(x, 0^, u^(x, t), t), 

6/j,(x,t) = Ty{ep^y{y,t) + ^/,^(t)); G M^^^ > 
0p,j,(y, t) = ^y(y, t)a^(y, t) - ^^.(y, t) (36) 

where 7^x(')5 "^yl') defined as in (fTHj) . and the functions \E'x(-)5 will specified later. 

Now we are ready to formulate the following result 

Theorem 2 (Properties of the interconnected systems) Let systems he given. 

Furthermore, suppose that the following conditions hold: 

1) The functions ipx{'^,t) , ■ipy{y,t) satisfy AssumptionU\ for systems / T^) . / T^) respectively; 

2) The systems 

Ipx = -^x{lpx, ^x, t) + C,x{t), Ipy = -'^yii'y, ^y, t) + Cy{t) (37) 

satisfy Assumption\^ with corresponding mappings 

xq; ^xj ||Cx(t)||2,[to,T]), 7y^,2(^l/o.^J/. IICj/(i)l|2,[to,T]), 

3) The systems have L\[tQ, oo] t-^ -^v2[to) oo] gains, that is 

||^x(x(t),t)||2,[to,T] < +7...2(IICx.(t)||2,[to,T]), 

||^,(y(t),t)||2,[to,T] < C^7. +7y,,.(IICy(t)l|2,[*„,T]), (38) 
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4) The functions fx{'x,6x,t), fy{y,6y,t) satisfy Assumptions\^^with corresponding con- 
stants Dx, -D^i; Dy, Dy^ and functions OLx{^.,t), OLy{y,t); 

5) The functions /;,^(x, y,t), hy{'x,'y,t) satisfy the following inequalities: 

||/i,(x,y,t)|| </3,||V^,(x,t)||, ||/i,(x,y,t)|| < /5,||^,(y,t)||, gM+ (39) 

Finally, let the functions \&^.(x, t), \&j^(y,t) in / T^) . / T^) satisfy Assumption\^ for systems 
|23); (Ell) respectively, and there exist functions pi{-), P2{-), Psi') > Id{-) G /Coo o'^c^ constant 
A e R+ s^xc/i i/ie following inequality holds: 

/^j/ ° 7^2.2° Pi ° + op30/3x07x.,.op20 (^-^ + 1^ (A) < A (40) 

for all A > A. Then 

CI ) The interconnection / t,'//'Jj) . / t,'/i)) with controls ^HB)> ^■^ forward- complete and trajec- 
tories x(t), y(t) are bounded 
Furthermore, 

C2) if properties hold for fx{'x,6x,t), fy{y,6y,t), /i^(x, y,t), /iy(x, y,t), anc? a/so 

functions (px{ilj^,LJx,t) , ipy{ipy,ujy,t), then 

lim V'x(x(t), t) = 0, lim V?y(y(i^), i) = (41) 

Moreover, 

Cc?^ if property holds for /^^(x, O^, t), fy{y, Oy, t), and the functions 
OL^ip^.t), dilJx{-K,t)/dt, OLyiy.t), diljyij ,t) / dt 
are locally hounded with respect to x, y uniformly in t, then 

lim /,(x(t), 61,, t) - /,(x(t), ^,(t), t) = 0, 
lim fyiyit), Oy, t) - fyiyit), ey{t),t) = 

r— >oo 

Let us briefly comment on the conditions and assumptions of Theorem El Conditions 1), 
2) specify restrictions on the goal functionals, similar to those of Theorem ^ Condition 3) 
is analogous to requirement to P3) in Theorem ^ condition 5) specifies uncertainties in the 
coupling functions hx{-)., hy{-) in terms of their growth rates w.r.t. ipxi'), '4'yi')- We observe here 
that this property is needed in order to characterize the L2 norms of functions /;,,.(x(t), y(t), t), 
hy{'x(t),'y(t),t) in terms of the L2 norms of functions ipx{'^(t),t), i'yiy(t),t)- Therefore, it is 
possible to replace requirement (jH^ with the following set of conditions: 

\\hxi^it),yit),t)h,[t„T] < /3x\\MAt),t)h,[t„T]+Cx, 
\\hy{^{t)Mt),t)h,lto,T] < PyUyiyit)M2,lto,T]+Cy 

The replacement will allow us to extend results of Theorem |21 to interconnections of systems 
where the coupling functions do not depend exphcitly on ipx{'^(t),t), ■ipy{y{t),t). We illustrate 
this possibility later with an example. 



11 



Condition pUj) is the small-gain condition with respect to the Ll[tQ,T] norms for intercon- 
nection (jnni), (pTTjl with control (jH^ . In the case that mappings 'yx2,2i')y 7»/2,2(') (|SZI) are 
majorated by linear functions 



7x2.2 (A) < fe,,A, 7,,_,(A) < gy, ,A, A > 0, 



condition (jlUj) reduces to the much simpler 

Mx9x2,29y2.2 + ^) + 1^ < 1 

Notice also that the mappings 7^2,2 (■)) 7y2,2(') are defined by properties of the target dynamics 
()37p. and, in principle, these can be made arbitrarily small. This eventually leads to the 
following conclusion: the smaller the L2-gains of the target dynamics of systems iSi, 1S2, the 
wider the class of nonlinearities (bounds for (3xi Py, domains of D^, -Di,a;, Dy, Di^y) which admit 
a solution to Problem ^ 

Example Let us illustrate application of Theorem |21 to the problem of decentralized control of 
two coupled oscillators with nonlinear damping. Consider the following interconnected systems: 

Xl =X2 f 2/1 = 2/2 .^^x 

X2 = fxiXl, 9x) + hyi + U^, \ 1/22 = fy{,yi, Oy) + ^23^1 +Uy, 

where fci, k2 eM. are uncertain parameters of coupling, functions f{xi,9x), f{yi,9y) stand for 
the nonlinear damping terms, and 6^, Oy are unknown parameters. For illustrative purpose we 
assume the following mathematical model for functions fx{-), fy{-) in ()44p : 

fx{xi, Ox) = Ox{xi - xq) + 0.5 sin(6'x.(xi - xq)), , ^ 

(45) 

fyiVi^ ^y) = ^y{yi " 2/o) + 0.6 sin(6'y(yi - ?/o)) 

where Xq, 2/0 are known. Let the control goal be to steer states x and y to the origin. Consider 
the following goal functions 

i)x{^, t)=Xi+ X2, ipyiy, t) =yi + y2 (46) 
Taking into account equations (jUj) and (jlUj) we can derive that 

xi = -xi + ?/'a;(x(t), t), yi = -yi + 'ipy{y, t) (47) 
This automatically implies that 

lki(t)||oo,[to,r] < \\xi{to)\\ + ||^^(x(t),t)||oo,[to,T] 

\\yi{t)\\oo,ito,T] < \\yi{to)\\ + \\^y{y{t),t)\\oo,[to,T] 

Hence, Assumption [T] is satisfied for chosen goal functions ipx{-) and ipy{-). Notice also that 
equalities (|17j) imply that 

\\xim2,[to,T] < 2-l/lxi(/o)|| + \\M^,t)h,[t„T] 
\\yi{t)h[to,T] < 2-'/'\\y,{to)\\ + ||V^,(y,t)||2,[*o,T] 
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Moreover, according to ()47|) limiting relations 



guarantee that 



lim ipx{'^it),t) = lim Xi{t) + X2{t) = 0, 
lim 7/>,(y(t), t) = lim + t/2(t) = 

t—>oo t— >oo 



lim Xi(t) = 0, lim X2(t) = 0, lim yi(t) = 0, lim y2(t) = 

t— >oo t— >oo t— >oo t— >oo 



Hence, property ()49|1 ensures asymptotic reaching of the control goal. 
According to equations (jHSl), (ISBI) control functions 

Ux = -Ki'x -X2- fx{Xl, Ox) 



(50) 

'^y = -^i^y - 2/2 - fyiVi, Oy), Xx, Xy>0 



transform system ()44|) into the following form 

i^x = -Xxi^x + fx{xi, 9x) - fx{xi, Ox) + hyi 

(51) 

= -Xxipx + /x(a;i, 6':r) - fx{xi, 9x) + ^2X1 

Notice that systems 

i^X = -Xxi^X + ^x{t), i)y = -Xytpt + ^yit) 

satisfy Assumption |21 with 

7x2,2 = ^||V'x(x(t),t)||2,[to,T], 7?/2,2 = ^\\i^y{yit),t)hlto,T] 

Ax Ay 

respectively, and functions fx{-), fy{-) satisfy Assumptions El IH with 

Dx = 1.5, Dx,i = 0.5, ax(x,t) = xi - xq, 
Dy = 1.6, Dy^i = 0.4, ay{y,t) = - 

Hence conditions l)-4) of Theorem |2l are satisfied. Furthermore, according to the remarks re- 
garding condition 5) of the theorem, requirements (jH^ can be replaced with implicit constraints 
fjl^ . These, however, according to PHj) also hold with Px = ^1, Py = ^2- 

Given that ^^^(x, t) = xi —Xq, ay{y, t) = yi — yo, Assumption El will be satisfied for functions 
ax{'x,t), ay{y,t) with \E'x(-) = 0, \E'y(-) = 0. Therefore, adaptation algorithms (jHUj) will 
have the following form: 

9x = Txiixi + X2)ixi - Xo) + 9xj), 

dx,I = Xx{xi + X2){xi - Xo) - {Xi + X2)X2 

9y = Ty{{yi+y2){yi-yo) + 9y,i), (52) 
9y,i = Xy{yi + 2/2) (2/1 - yo) - (2/1 + 2/2)2/2 

Hence, according to Theorem El boundedness of the solutions in the closed loop system (jSIJ, 
(jK^ is ensured upon the following condition 
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Figure 2: Plots of trajectories Xi{t) (panel a), X2(t) (panel b), yiit) (panel c), y2it) (panel d) as 
functions of t in closed loop system (jH)), (jSUjl . (jH2I)- Dotted lines correspond to the case when 
ki = k2 = 0.4, and solid lines stand for solutions obtained with the following values of coupling 
= 1, ^2 = 0.1 



Moreover, given that properties 1^- HSlhold for the chosen functions '?/'^(x, t), il)y{y,t), condi- 
tion (j53|l guarantees that limiting relations (j4T| . (jl2|l hold. 

Trajectories of the closed loop system ()44j). (j^Uj) . ()52|) with the following values of parameters 
Vx = Ty = 1, \x = \ = 2, a;o = Z/o = ^'j,. = 6*2^ = 1 and initial conditions Xi(0) = — 1, 
X2(0) = 0, 2/i(0) = 1, 2/2(0) = 0, 4,/(0) = -1, ^y,/(0) = -2 are provided in Fig. H 

5 Conclusion 

We provided new tools for the design and analysis of adaptive decentralized control schemes. 
Our method allows the desired dynamics to be Lyapunov unstable and the parametrization of 
the uncertainties to be nonlinear. The results are based on a formulation of the problem for 
adaptive control as a problem of regulation in functional spaces (in particular, L\[tQ, T] spaces) 
rather than of simply reaching of the control goal in R"^. This allows us to introduce adaptation 
algorithms with new properties and apply a small-gain argument to establish applicability of 
these schemes to the problem of decentralized control. 

In order to avoid unnecessary complications, state feedback was assumed in the main-loop 
controllers which transform original equation into the error coupled model. Extension of the 
results to output-feedback main loop controllers is a topic for future study. 
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6 Proofs of the theorems 



6.1 Proof of Theorem [T] 

Let us first show that property PI) holds. Consider solutions of system (jH), (fTUj) . (fTTj) . (fTHj) 
passing through the point x(to), Gi{to) for t G [to;^*] • Let us calculate the time-derivative of 
function ^(x, t): 0{x, t) = T{Op + Oj) = T{i)oc{x, t) + ?/'q;(x, t) - ^(x, t) + Oj). Notice that 

^a(x, t) - *(x, t) + 61, = V^(x, t) ^ ^ xi + 7/;(x, t)— ^X2+ 

dXi dX2 .g^x 

^, ga(x,t) d^{^,t) . gvl/(x,t) . cM(x^ A ^ ^ 

^^^^ 9^''^ " " + ' 

According to Assumption ^'^J^f^ = '?/'(x, t) • Then taking into account (jSH), we obtain 

9xi 9X1 



^d(x, t) - ^(x, t) + 6>, = ( ^(x, t)^^^^ - 1^ ) xi 



(55) 



dt dt 

Notice that according to the proposed notation we can rewrite the term ^^/'(x, t)^^^ — ^ j xi 
in the following form: ^/^(x, )f:)Lf^Q;(x, t) — Lf^\E'(x, t) + ('?/'(x, )f:)Lg^a(x, t) — Lg^\l'(x, t)) m(x, 0, t). 
Hence, it follows from p7|l and that tljct{'x,t) — \E'(x, t) + 0j = (p{tp)a{'x,t). Therefore, the 
derivative 6{'x,t) can be written in the following way: 

6/ = r(^ + (/;(^))a(x,t) (56) 

Asymptotic properties of nonlinear parameterized control systems with adaptation algorithm 
fl5(i|l under assumption of Lyapunov stability of the target dynamics were investigated in j22|- In 
the present contribution we aim to provide characterizations of the closed loop system in terms 
of functional mappings between functions '?/'(x(t), t), e(t), and /(x(t), 6, t) — /(x(t), 0(t),t) and 
without requiring Lyapunov stability of the target dynamics flTT|l . 
For this purpose consider the following positive-definite function: 



v^{e,9,t) = -\\e-er^-, + ^j^ s\T)dT (57) 

Its time-derivative according to equations can be obtained as follows: 

v^,{e,e,t) = {v>W + i^){e - 0)^«(x,t) - ^^e\t) (58) 

Hence using Assumptions El ID and equality (jlOj) we can estimate the derivative Vg as follows: 

v^{e, e, t) < -(/(x, e, t) - /(x, e, t) + eme - efa{^, t) - -^^e\t) 
< -;^(/(x, e, t) - /(x, e, t)f + I l/(x, e, t) - /(x, e, t) \ 

^ie\t) < -1 f|/(x,0,t) -/(x,0,t)| - ^e{t)\ < (59) 
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It follows immediately from (jFfj) that 

\\e{t)-o\\l., < \\e{to)-d\ 



2 4- 

r-i + 



D 

2Dl 



2 

2,[to,oo] 



(60) 



^11^1 + 



In particular, for t G [to?^*] we can derive from ()57p that \\0{t) — ^||p-i < ||^(^o' 

. Therefore 9{t) G i:^[to,T*]. Furthermore |/(x(t), 6>(t), t) - /(x(t), 0, t)| 



2§^£:(t) G L2[to,7'*]. In particular 



|/(x(t),^(t),t)-/(x(t),0,t)| 



< 



2,[to,T*] 



^-^(to)llr-i + 



4Df 



2 



(61) 



Hence /(x(t), ^(t), t) - /(x(t),0,t) G L^[to,T*] as a sum of two functions from L^[to,T*]. In 
order to estimate the upper bound of the norm ||/(x(t), ^(t), t) — f{-x(t),6,t)\\2^[to,T*] from (jHH) 
we use the Minkowski inequality: 

D 



/(x(t),0(t),t)-/(x(t),0,t)| 



2L)i 



< 



^\\O-0ito) 



0.5 



+ 



D 

2lTi 



and then apply the triangle inequality to the functions from Ll[to,T* 



||/(x(t),^(t),t)-/(x(t),0,t)||2,[*o,T1< 

/(x(t),^(t),t)-/(x(t),0,t)--^£(t) 



2Di 



+ 



(62) 



D 

2lTi 



2,[to,T*: 



< 



D 



\e-e{t. 



0.5 



2,[to,T*] 

D_ 



2,[to,T*] 



Therefore, property PI) is proven. 

Let us prove property P2). In order to do this we have to check first if the solutions of 
the closed loop system are defined for all t G IR+, i.e. they do not go to infinity in finite 
time. We prove this by a contradiction argument. Indeed, let there exists time instant tg such 
that ||x(t,)|| = oo. It follows from PI), however, that /(x(t), ^(t), t) - /(x(t),0,t) G Ll[to,ts]. 
Furthermore, according to (jU^ the norm ||/(x(t), 6(t),t) — /(x(t), ^, t)||2,[to,ts] can be bounded 
from above by a continuous function of 6, ^(to)? T, and ||£^(t)||2,[to,oo]- Let us denote this 
bound by symbol Df. Notice that Df does not depend on tg. Consider system (fTUI) for 
t G [to,ts]: ^ = /(x,6',t) - /(x,^,t) - ip{ip,LJ,t) + e{t). Given that both /(x(t),0,t) - 
/(x(t), 6{t),t),e{t) G Ll[tQ,ts] and taking into account Assumption |21 we automatically obtain 
that ip{'x(t),t) G Ll^[tQ,ts]. In particular, using the triangle inequality and the fact that the 
function 700,2 (V'(xo, ^o); ^5 ^) in Assumption |21 is non-decreasing in M, we can estimate the 
norm ||^(x(t), t) ||oo,[to,t,] as follows: 



||^(x(t),t)||oo,[to,t.] < 7oo,2 {ipiy^o,to),u),Df 



2, [to, 00 



) 



(63) 
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According to Assumption ^ the following inequality holds: 

||x(t)||oo,[to,t.] < 7 (Xo,6',7oo,2 (^(Xo,to),^,^/ + Mt)\\l[to,oo])) (64) 

Given that a superposition of locally bounded functions is locally bounded, we conclude that 
||x(t)||oo[to,is] is bounded. This, however, contradicts to the previous claim that ||x(ts)|| = oo. 
Taking into account inequality (jHIH) we can derive that both 0{'x{t),t) and Oj{t) are bounded 
for every t G IR+. Moreover, according to (jMjl . (jHUIl these bounds are themselves lo- 

cally bounded functions of initial conditions and parameters. Therefore, x(t) G L^[to, oo], 
0{'x{t),t) G L'^[to,oo]. Inequality (PT|) follows immediately from (jHH), (US), and the triangle 
inequality. Property P2) is proven. 

Let us show that P3) holds. It is assumed that system (fTT|) has Ll[tQ,oo] t-^ Lp[tQ,oo], 
p > 1 gain. In addition, we have just shown that /(x(t), 6, t) — f{-x(t), 0(t),t),e(t) G L2[to, oo]. 
Hence, taking into account equation (fTUjl we conclude that ip{'x{t),t) G Lp[to, oo], p > 1. On the 
other hand, given that f{ii.,d,t), ip{ip,uj,t) are locally bounded with respect to their first two 
arguments uniformly in t, and that x(t) G LJ^[to, oo],'?/'(x(t), t) G L^[to,oo], 0{t) G LJ^[to,oo], 
d G fie, the signal (^(^(x(t), t), cu, t) + /(x(t),0,t) - /(x(t), t) is bounded. Then e{t) G 
L^[to, oo] implies that ip is bounded, and P3) is guaranteed by Barbalat's lemma. 

To complete the proof of the theorem (property P4) consider the time-derivative of function 
/(x,^,t): 

^/(X,^,t) = ^f(x,0)+g(xMx,0,t)/(X' ^>^) + 

^ ^{fW, ^, t) + V^)a(x, t) + — 

Taking into account that the function /(x, 6, t) is continuously different iable in x, 0; the deriva- 
tive df{'x,6,t)/dt is locally bounded with respect to x, uniformly in t; functions q:(x, t), 
dijj{'x, t)/dt are locally bounded with respect to x uniformly in t, then d/dt{f{'x, 0, t) — /(x, 0, t)) 
is bounded. Then given that f{'x(t),0,t) — f{'x(t),0(t),t) G Ll[to,oo] by applying Barbalat's 
lemma we conclude that /(x, 0, r) — /(x, ^, r) ^ as t ^ oo. The theorem is proven. 

6.2 Proof of Theorem [2] 

Let us denote 

A/,.[to,T] = ||/.(x,6>,,t)-/,(x,6l,,t)||2,[to,T], 

Afy[to,T] = \\Uy,0y,t)-fy{y,0y,t)\\2,[to,T]. 

As follows from Theorem^ the following inequalities hold 

Af4to,T]<C, + ^\\hy{^{t),y{t),t)h,ito,T] (65) 

^fy[to,T] < Cy+^\\h.i^it)Mt),t)\\2,ito,T], (66) 
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where C^, Cy are some constants, independent of T. Taking estimates (jUBj) into account 
we obtain the following estimates: 



A/.[to,T] + ||/i,(x(t),y(t),t)||2,[io,T] < 



c.+ i^ + i] \\hyi^{t)Mt)M2,[to,T] ^^^^ 



Afy[to,T] + ||/i,(x(t),y(t),t)||2,[to,T] < 

Dy , ,\ „, . „„ (68) 



The proof of the theorem would be complete if we show that the L2[to, T] norms of /i^(x(t), y(t), t), 
hy{':si{t),'y{t),t) are globally bounded uniformly in T. Let us show that this is indeed the case. 
Using the widely known generalized triangular inequality ^U] 

7(a + b)< 7((p + Id){a)) + 7((p + Id) o p-\b)), a, 6 G R+, 7, P e /Coo, 

equations (|^. (jUH|) and also property (jH^ . we conclude that 

||/i,(x(t),y(t),t)||2,[i„,r] < 

/3y ■ 7^2,2 ° Pi ((^ + 1) ||/ix(x(t),y(t),t)||2,[to,T]^ +C^,i 

||/l,(x(t),y(t),t)||2,[to,T] < 

Px- 1x2,2 °P2 ((;^ + ^) II^J/(x(t),y(t),t)||2,[to,T]^ +C'^.,1 



(69) 



where pi{-), P2(') ^ /Coo, Pi(')5P2(') > Id{-). Then, according to (jSHI), the existence of p3(-) G 
/Coo > -^c^(")5 satisfying inequality 

Py ° 7,2.2 ° Pi ° + 1^ o p3 o o 7,^^^ o p2 o (^-^ + 1^ (A) < A V A > A 

for some A G M+ ensures that the norms 

\\hy{^{t),y{t),t) \\2,[to,T], ||/ix(x(t), y(t), t) ||2,[io,T] 

are globally uniformly bounded in T. The rest of the proof follows from Theorem ^ The 
theorem is proven. 
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